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Abstract 

We analyse the loop induced production of Higgs boson pairs at future high-energy 
e~^e~ colliders, both in the Standard Model and in its minimal supersymmetric 
extension. The cross sections for Standard Model Higgs pair production through 
W/Z boson loops, e"'"e~ —>■ H^H^, are rather small but the process could be visible 
for high enough luminosities, especially if longitudinal polarization is made avail- 
able. In the Minimal Supersymmetric Standard Model, the corresponding processes 
of CP-even or CP-odd Higgs boson pair production, e~^e~ hh, HH, Hh and 
e~^e~ — > AA have smaller cross sections, in general. The additional contributions 
from chargino/neutralino and slepton loops are at the level of a few percent in most 
of the supersymmetric parameter space. 



*Supported by Deutsche Forschungsgemeinschaft DFG (Bonn). 



1. Introduction 



The search for scalar Higgs particles and the exploration of the electroweak gauge sym- 
metry breaking sector will be one of the main goals of future high-energy colliders. While 
Higgs particles will probably be first produced at the Large Hadron Collider LHC, the 
clean environment of high-energy e^e~ linear colliders would be required for a detailed 
investigation and for the high-precision measurement of the fundamental properties of 
the Higgs particles [Q, a crucial requirement to establish the Higgs mechanism as the 
basic mechanism to generate the masses of the known particles. To this end, a precise 
calculation of the branching ratios of all important decay channels as well as the cross 
sections of the production mechanisms, is mandatory 

The main production mechanisms of the Standard Model (SM) Higgs particle 
at high-energy e^e~ colliders are the bremsstrahlung process 0, e^e~ H^Z, 
and the vector boson fusion mechanisms e~^e~ — >■ W*W* / Z* Z* H^vv / H^e^ e~ . 
Besides providing the experimental signals for the Higgs particles, these processes would 
also allow to measure the Higgs couplings to gauge bosons, and to test that they are 
indeed proportional to the W/Z masses, a fundamental prediction of the Higgs mechanism. 
Some of the couplings of Higgs particles can also be measured by considering the decay 
branching ratios or higher-order production processes. For instance, the Higgs boson 
self-coupling can be accessible in the pair production processes e^e~ H^H^Z and/or 
W*W* /Z*Z* H^H^ 0], while the Yukawa coupling of light Higgs bosons to top quarks 
can be directly measured in the process e^e" —>■ tiH^ [0]. 

Higher order processes might also be useful to discriminate between the Higgs sector 
of the Standard Model from the usually more complicated scalar sectors of its possible 
extensions. Among these extensions, supersymmetric theories and in particular the Min- 
imal Supersymmetric Standard Model (MSSM) are the most natural, theoretically. In 
the MSSM, the Higgs sector is enlarged to contain two scalar doublet fields leading to a 
spectrum of five Higgs particles: two neutral CP-even {h and H), a CP-odd {A) and two 
charged (-f^^) Higgs bosons [0. In the decoupling regime |Q, where the A,H and if^ 
bosons are very heavy, the lightest Higgs boson h has exactly the same properties as the 
SM Higgs boson, except that its mass is restricted to be smaller than Mh < 140 GeV. 

If the genuine supersymmetric particles were too heavy to be kinematically accessi- 
ble in collider experiments, the only way to distinguish between the SM and the lightest 
MSSM Higgs boson in this limit, is to search for loop induced contributions of the su- 
persymmetric particles, which could give rise to sizeable deviations from the predictions 
of the SM. Well known examples of this loop induced processes are the 77 widths of the 
Higgs particles [which can be measured in the decay Higgs —>■ 77 at hadron colliders, and 
more precisely at high-energy e~^e~ colliders in the direct production of Higgs particles via 
laser-photon fusion, 77 Higgs] or the process e'^e~ Z+Higgs which in the MSSM 
receive extra contributions from supersymmetric gaugino and sfermion loops 0, |ll[ . 



Another type of such discriminating processes is the pair production of Higgs bosons 
which will be analyzed in this paper. In the Standard Model, where it has been first 
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discussed in Ref.[|T2|, the process 

e+e- ifOi/O 

is mediated only by W and Z boson loops, Fig. la, while in the Minimal Supersymmetric 
extension, additional contributions to the corresponding process 

e+e" hh (2) 

will originate from chargino, neutralino, selectron and sneutrino loops, as well as loops 
built up by the associated A and bosons; Fig. lb. At high center-of-mass energies, 
^/s > 1 TeV, the cross sections of the processes are rather small, being of the order of a 
fraction of a femtobarn. However, with the large luminosities expected at these colliders, 
J C > 100 fb~^, and with the availability of longitudinal polarization of the electron 
and positron beams which could increase the cross sections by a factor of four, a few 
hundred events might eventually be collected in the course of a few years, allowing for 
the experimental study of this final stateQ. 

The supersymmetric contributions to the process e~^e~ — > hh in the MSSM turn out 
to be rather small compared to the dominant W/Z [and A/H^] contributions. Except in 
some regions of the SUSY parameter space where they can reach the level of ~ 15%, they 
are typically of the order of a few percent, rendering the distinction between the SM and 
the light MSSM Higgs bosons rather difficult in the decoupling regime. For completeness, 
we have also analyzed the pair production of the heavy neutral Higgs particles of the 
MSSM. Up to phase space suppression factors, the cross sections for the processes 

e+e" ^ HH and AA (3) 

are of the same order as the SM Higgs pair production cross section. In contrast, due to 
mixing angle factors suppression, the cross section for 

e+e- ^ Hh (4) 

is at least one order of magnitude smaller, but the supersymmetric contributions are 
of the same size as the W/Z /Higgs contributions. The processes e+e~ — » HA/hA and 
e+e~ H^H'^ appear already at the tree level, and will not be considered hereQ. 

The rest of the paper is organised as follows. In the next section we present the cross 
sections for the Higgs pair production process in the Standard Model. In section 3, we 
discuss the case of the Higgs pair production process in the MSSM. Some conclusions will 
be drawn in section 4. In the Appendix, we will give the lengthy analytical expressions 
of the cross sections in the MSSM. 



^This process is similar to the double Higgs production W*W* H'^H^ or e+e^ ZH^H^ 
vvH'^H^ , which consist of the same final state [although with missing energy] and which can be used to 
measure the trilinear Higgs self-coupling 

^The radiative corrections to these processes, including the contributions of the supersymmetric par- 
ticles, have been recently discussed in Refs. 10, 13 . 
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2. SM Higgs Pair Production 

In the Standard Model, many Feynman diagrams could in principle contribute to the 
amplitude for the Higgs boson pair production process, eq.(l). However, in the limit of 
vanishing electron mass, which implies a vanishing Hee coupling, the contributions of 
the tree-level diagrams where the Higgs bosons are emitted from the electron lines are 
negligible. Non-zero contributions to the process e^e~ —>■ H^H^ can therefore only come 
from one-loop diagrams. 

Because of the exact chiral symmetry for = 0, all diagrams involving the one-loop 
ife^e^ vertex must give zero contributions [in fact this statement is valid to all orders 
in perturbation theory]. Furthermore, because of CP invariance, the amplitudes of the 
diagrams with 7 and Z boson s-channel exchanges, which give rise to two Higgs bosons 
through vertex diagrams, also vanish [the contribution of the longitudinal component 
of the Z boson will be proportional to the electron mass, and is therefore negligible]. 
Additional contributions from vertex diagrams involving the quartic WW / Z Z 
couplings are proportional to the electron/neutrino masses and are again negligible. The 
only contribution to Higgs pair production in the Standard Model will therefore come 
from W and Z box diagrams. Fig. la. 

Working in the Feynman gauge, where the W/Z boson contributions can be split 
into the g^y and the neutral Goldstone parts, the ultraviolet-finite amplitudes have been 
reduced from a complicated tensorial form to scalar Passarino-Veltman functions 



using the package FORM |T3[. Before the reduction to scalar functions, we have verified 



that our expressions agree with those obtained in Ref. ||12||. The amplitudes have then 
been symmetrized and squared; a factor of 1/2 has been included because of the two 
identical particles in the final state. Allowing for the polarization of both the electron 
and positron beams, which are respectively denoted by A_ and A+, the differential cross 
section of the process (1) [in the limit me = 0] is then given by 



dcos^ 1024 7r5 s{ut - Mfj] 



(1 + A_)(1-A+) 



pvv _|_ \- =W ^) pZ 



+ (l-A_)(l + A^)4jF^n (5) 

with s, t and u the usual Mandelstam variables, 9 the angle between the initial electron 
and one of the Higgs bosons, and s^r = 1 — = siv? 6w The form factor with 
V = W/Z can be split into two terms corresponding to the contributions of the (yf^j, and 
Goldstone parts 

F^ = 2M^F^ + F^ (6) 

where, in terms of the scalar two, three and four-point Passarino-Veltman functions^, 
■^The complete analytical expressions of the scalar functions can be found in Ref.pq] for instance; for 



their numerical evaluation we have used the package FF |17 
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denoted respectively by Bq, Cq and Dq, the form factors F^^ given by 



V 



with 



2M]jU-My{u-t)-u{u + t) 
+ {Ml - u) [Co , 0, , 0, u, Ml 



- 1) Co {Ml 



M^, 0, 0, s) 
j)+Co{0,MlM^,0,Ml 



tl,Ml,s,Ml,Ml) \ + \t 



* — > u , 



(7) 



I [tu - + 2u Ml - 2 Mfj] Co (m^, 0, M^, 0, 0, s) 

+u (u - Ml) [Co (m^, 0, M^, 0, u, Ml) + Co (o, M^, M^, 0, Ml, u) 
t^u-v?t-2u^ + 4.Ml(tu-v?) + Miit + u) ^ /, . , . , . 

+ [1*^ - + + -u) + 2?/M| {Ml -u)-2Ml M^Dq 

+f^^^[Bo{M'y,M'y,Ml)-Bo{M^y,M'y,s)\^ + |^ — 



t) + Ml,u(t + u) + 2uMl{Ml - u) - 2MlM^ 
\)-Bo (m2,M^,s) 



^0 



Do = Do{Ml, 0, M^, Ml, 0, 0, Mf, M^, s, u) 



(8) 



(9) 



The cross sections are shown in Fig. 2a as a function of the Higgs boson mass for two 
center-of-mass energies, ^/s = 500 GeV and 1.5 TeV. Except when approaching the 2Mh 
threshold [and the small dip near the WW threshold], the cross sections are practically 
constant for a given value of the cm. energy, and amount to cr ~ 0.2 fb at ^/s — 500 GeV 
in the unpolarized case. The decrease of the cross sections with increasing center-of-mass 
energy is very mild: at ^/s = 1.5 TeV, the cross section is still at the level of cr ~ 0.15 fb 
for Higgs boson masses less than Mh ^ 350 GeV. 

The effect of polarizing the initial electron/positron beams is also shown in Fig. 2a. 
With left-handed polarized electrons, the cross section e^e"*" H^H^ is larger by a 
factor of two, while for left-handed electrons and right-handed positrons, the cross section 
e^e^ H^H^ is larger by a factor of four, compared to the unpolarized case. [These 
simple factors of 2 and 4 are due to the fact that the contribution of the W boson is much 
more important than the one of the Z boson, a mere consequence of the larger charged 
current couplings compared to the neutral current couplings.] Therefore, the availability 
of longitudinal polarization of the initial beams is very important. 

With integrated luminosities of the order of / £ ~ 100 fb^^ which are expected to be 
available for future high-energy linear colliders, one could expect a few hundred events in 
the course of a few years, if both initial beams can be longitudinally polarized. 
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Fig. 2b shows the angular distributions dcr/dcos^ at a center-of-mass energy of = 
500 GeV for a Higgs boson mass of Mh = 150 GeV. It is forward-backward symmetric, 
a consequence of the two identical particles in the final state, and follows approximately 
the dcr/dcos6' ~ sin^ 6* law. Again, the angular distribution does not strongly depend on 
the Higgs mass and on the cm. energy, except near the 2Mh production threshold. 

For Higgs masses below Mh ^ 140 GeV, the signal will mainly consist of four b quarks 
in the final state, 

e+e- ^ H^H^ hi hi (10) 

since the dominant decay mode of the SM Higgs boson in this mass range is — > hh. This 
calls for very efficient /x-vertex detectors to tag the h quark jets. Since these rare events 
will be searched for only after the discovery of the Higgs boson in the main production 
processes, the Higgs boson mass will be precisely known and the two mass constraints 
m{hl) = Mh, together with the rather large number of h quarks in the final state, give a 
reasonable hope to experimentally isolate the signals despite of the low rates. 

For larger Higgs masses, Mh > 140 GeV, since W'^W" and H° ZZ will 

be the dominant decay modes of the Higgs boson, the signals will consist of four gauge 
bosons in the final state, 

e+e- ^ ^ W+WW^W- , W^W'ZZ , ZZZZ (11) 

leading to eight final fermions. These rather spectacular events should also help to exper- 
imentally isolate the signal. 

Although the discussion of the background events is beyond the scope of this paper, 
a few comments are nevertheless in order. 

{i) For light Higgs bosons, assuming that b quarks will be efficiently tagged, the main 
backgrounds will consist of the QGD 4 6-jet process, ZH'^ and ZZ pair production with 
the Z bosons decaying into hi pairs, the latter background being the most dangerous 
especially if Mh ~ Mz- However, since ZZ production is mediated by t-channel electron 
exchange, the cross section is strongly peaked in the forward and backward directions 
contrary to the signal process as shown in Fig. 2b; a strong cut on cos 6 should considerably 
reduce the background, while leaving the signal almost unaltered. 

(ii) For higher Higgs boson masses, multiple vector boson production will be the main 
background; however, since these are higher-order processes in the electroweak coupling, 
the cross sections should be small enough for this background to be manageable. 

{Hi) Finally, one should note that double Higgs production in 1^1^ fusion, e'^e~ — > 
H^H^vv and in the bremsstrahlung process e+e~ — > H^H^ + Z[-^ uP] would also act as 
backgrounds since the two neutrinos will escape undetected. However, the requirement 
of no missing energy in the final state will easily discard these events. 
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3. MSSM Higgs Pair Production 

In the case of the pair production of the CP-even Higgs bosons of the Minimal Supersym- 
metric Standard Model, e+e^ $i$2 with $i,$2 = h,H, several additional Feynman 
diagrams will contribute to the processes; Fig. lb. Besides the W and Z boson box di- 
agrams, one first has the box diagrams with the exchange of the pseudoscalar and the 
charged Higgs bosons, A and H^. Then, one has two classes of box diagrams built up by 
chargino/sneutrino and neutralino/selectron loops^. These two sets consist of diagrams 
where both Higgs bosons couple to the neutralino/chargino or to slepton pairs, and dia- 
grams where one Higgs boson couples to neutralinos/charginos and the other to sleptons. 
One has also to include the crossed diagrams in Fig. lb, with the exchange of $i and $2- 

For the pair production of the CP-odd Higgs boson in the MSSM, e~^e~ AA, only a 
small subset of the previous Feynman diagrams contributes; Fig.lc and the corresponding 
crossed diagrams. Indeed, because of CP invariance the pseudoscalar Higgs boson does 
not couple to vector boson pairs and to slepton pairs [in the latter case, the couplings 
would be proportional to the lepton masses, rrie and m^]. The only relevant contributions 
will therefore come from the diagrams with h, H, boson exchange, and those with 
chargino/sneutrino and neutralino/selectron exchanges, where both A bosons are emitted 
form the neutralino/chargino lines. 

We have calculated the cross sections of the four processes eqs.(2,3) and (4) in the 
MSSM, following the same steps as those discussed in the case of the SM Higgs pair 
production. The analytical expressions of the cross sections are much more involved than 
in the SM case, a consequence of the many additional contributions, and in the case 
of e~^e~ — > Hh, of the two different particles in the final state. These expressions will 
therefore be given in the Appendix. Here, we will simply describe our inputs and discuss 
the numerical results which have been obtained. 

Besides the four masses Mh, M^, Ma and Mjj± of the Higgs particles, two additional 
parameters determine the Higgs sector of the MSSM at the tree level: the ratio tan j3 = 
V2/V1 of the vacuum expectation values of the two Higgs doublet fields and a mixing 
angle a in the neutral CP-even sector. However, supersymmetry leads to several relations 
among these parameters and only two of them are in fact independent : if the pseudoscalar 
Higgs mass and tan/5 are specified, all other masses and the mixing angle a can be 
derived at the tree-level [jl|. Supersymmetry imposes a strong hierarchical structure on 
the mass spectrum, Mh < Ma < Mh, Mw < Mh± and Mh < Mz, which however is 
broken by radiative corrections 



The leading part of the radiative corrections grows as the fourth power of the top 



quark mass rrit and the logarithm of the common squark mass Ms [|18] . At the subleading 



level the radiative corrections will introduce the supersymmetric Higgs mass parameter /i 



and the parameters At, Aj, in the soft symmetry breaking interaction |T9[. In our analysis 



'^Vertex diagrams involving the quartic slepton-slepton-$i-'I>2 couplings will give rise to contribu- 
tions which are proportional to the mixing between the two slepton eigenstates. Since this mixing is 
proportional to the partner lepton masses rrig and m^, these diagrams give negligible contributions. 
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we will take into account the full one-loop corrections for the Higgs boson masses ; for 
the mixing angle we will use the one-loop improved pO| relation [tan/5 and are the 



input parameters and is the full one-loop corrected h boson mass] 

-(M^+M|) tan/? , , 

tana = ^-A^ ^ -^r- 12 

M| + Ml tan2 /3 - (l + tan2/3)Mft ^ ' 

Once these parameters are fixed, the couphngs of the Higgs particles to gauge bosons 
and fermions will also be uniquely fixed; the Higgs-vector boson couplings relevant to our 
analysis are given in Table la. 

To fully describe the supersymmetric sector two parameters must be introduced, in 
addition to M^, tan/5, M5, At, Af, and /i. These are the gaugino mass parameter M2 
[we will use the GUT relation which fixes the bino mass in terms of the gaugino mass. 
Ml = I tan^ 6w M2] and the common slepton mass M^. The masses of the four neutralinos 
Xi [i = 1, ••,4] and the two charginos xt = l^^] are then completely fixed by tan/?, 
and M2. The masses of the left- and right-handed selectron e^,, e/j, and of the left-handed 
electronic sneutrino 9^ will be fixed by Mi [in practice these masses will approximately 
be given by mg^ ~ rrig^ ~ m£>^ = M^; the mixing between the two selectron states is 
negligible] . 

For the couplings of the Higgs bosons to neutralinos, charginos and sleptons, one also 
will need the value of the mixing angle a which is fixed by tan (3 and Ma- These couplings 
are given in Tables Ib-lc; the couplings of the electrons to charginos/neutralinos and 
sleptons which will also be needed in the analysis are given in the Appendix. 

We are now in a position to discuss the numerical results. The cross sections for the 
four processes e~^e~—>-hh, HH, AA and hH are displayed in Figs. 3-6, as a function of the 
Higgs boson masses and for two values of the cm. energy ^/s = 500 GeV and 1.5 TeV 
[except for hH production] and two representative values of tan (3 = 1.5 and 50. In all the 
figures we have chosen the parameters M2 = —fi = 150 GeV, while the common slepton 
and squark masses are taken to be Ml = 300 GeV and Ms = 500 GeV; the parameter At 
and At) are set to zero. Only the unpolarized cross sections are discussed: as mentioned 
previously, they are simply increased by a factor of 2(4) when the initial beams(s) are 
longitudinally polarized. 

Fig. 3 shows the cross section for the process e'^e^—>-hh. The solid lines are for the full 
cross sections, while the dashed lines are for the cross sections without the SUSY contri- 
butions [this will correspond to a two-Higgs doublet model with the MSSM constraints]. 
Let us first discuss the case where the supersymmetric contributions are not included. 
For small tan f3, the cross section is of the same order as the SM cross section and does 
not strongly depend on M/j especially at very high-energies. Although the WWh/Zhh 
couplings are suppressed by sin(/5 — a) factors, the suppression is not very strong and the 
W/Z box contributions are not much smaller than in the SM; the diagrams where A/H^ 
are exchanged will give compensating contributions since the hAZ/hH^W couplings are 
proportional to the complementary factor cos(/3 — a). As in the SM case, the cross sections 
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slightly decrease with increasing energy. 

For large tan /3 values, the factors sin / cos(/3 — a) vary widely when is varied. For 
small Mfi, the factor sin(/9 — a) — > 0, and the contribution of the diagrams with A/H^ 
exchange dominates. The latter contribution decreases with increasing Al^ [i.e. with 
decreasing cos(/9 — a)], until the decoupling limit is reached for Mh — 110 GeV. In this 
case, the factor sin(/3 — a) ^ 1 and the W/Z boson loops are not suppressed anymore; 
one then obtains the SM cross section. 

The contributions of the chargino/selectron and neutralino/sneutrino loops lead to a 
destructive interference. At high-energies, the supersymmetric boxes practically do not 
contribute; but at low energies, and especially below the decoupling limit, the SUSY 
contributions can be of the order of ~ 10%. We have scanned the SUSY parameter 
space, and the maximum contribution of the SUSY loops that we have found was about 
~ —15%. In the decoupling hmit, the SUSY contributions are, at most, of the order of 
a few percent. Because of the rather low production rates, it will therefore be difficult to 
experimentally see this effect. 

Fig.4 and 5 show the cross sections for the processes e+e~ — > HH and e+e~ — > AA, 
respectively. Except for small A/H masses where they are of the same order as in the 
SM, the cross sections are below the 0.1 fb level and the signals will be hard to detect, 
especially for Ma-iMh ^ 350 GeV [independently of the phase space suppression]. The 
SUSY contributions are relatively even smaller than for the case of e+e~ hh. 

Finally, Fig. 6 shows the cross section for the case e~^e^ — > hH at y/s = 500 GeV. It is 
one order of magnitude smaller than in the previous cases and therefore negligibly small. 
This is due to the fact that the W/Z and A/H^ contributions are both suppressed since 
the combination sin(/? — a) x cos(/? — a) appears in all these contributions; since one of 
the factors is always small, this brings these contributions down to the level of the small 
SUSY-loop contributions. 

4. Summary 

We have analyzed the one-loop induced production of Higgs boson pairs at future high- 
energy e+e^ colliders in the Standard Model and in the Minimal Supersymmetric Standard 
Model. In the SM, the unpolarized cross section is rather small, of the order 0.1-0.2 fb. 
The longitudinal polarization of both the e~ and e"*" beams will increase the cross section 
by a factor of 4. With integrated luminosities J C > 100 fb^^ as expected to be the case 
for future high-energy linear colliders, one coTild expect a few hundred events in the course 
of a few years if longitudinal polarization is available. The final states are rather clean, 
giving a reasonable hope to isolate the signals experimentally. 

In the MSSM, additional contributions to the processes e+e" hh, HH, Hh and 
e'^e" — > AA come from chargino/neutralino and slepton loops. For hh production, the 
contributions of the supersymmetric loops are in general rather small, being of the order 
of a few percent; the cross sections are therefore of the same order as in the SM. For 
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the processes involving heavy Higgs bosons, the cross sections are even smaller than for 
e'^e" — > hh, and the signals will be hard to be detected experimentally. 



Acknowledgements: We thank W. HoUik for discussions and J. Rosiek for providing 
us with the Fortran code for the one-loop corrected MSSM Higgs boson masses. 



APPENDIX: Cross sections in the MSSM 

In this Appendix, wc will give the analytical expressions of the cross sections for the pair 
production of the MSSM Higgs bosons. We start with the production of two CP-even 
Higgs bosons, e'^e" — > $i$2 with $i, $2 = h, H. Similarly to the SM case, the differential 
cross section for the process can be written as: 



da 



1 



w 



dcos^ 1 + (5$i$2 Sn^s 
(1 + A_)(1-Ah 



2 \2 
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k=l 



1/2 



{ut-M^Ml] 



+ (1-A_)(1 + A+) 



fc=i 



(Al) 



where A_ and A+ are the longitudinal polarizations of the initial e~ and beams, 9 the 
scattering angle and s,t,u the usual Mandelstam variables; 5(i>j#2 = 1 in the case where 
we have two identical Higgs bosons in the final state, otherwise 5$;^$2 = 0. 

For the class of diagrams consisting of W/Z and A/ loops [a two Higgs-doublet 
model], there are three contributing amplitudes Fi 2,z given by 



Fa{Mz, 0, Mz, Mz, s, t) + Fa{Mw, 0, Mw, Mw, s, t) 



w 



Fr = g^,vv9^,vv'-^r^ FAiMz,0,Mz,Mz,s,t) 



--w 



9^\VV 9^ivv 



(2s 



Fb{Mz, 0, Mz, Mz, s,t) + \ Fb{Mw, 0, M^, M^, s, t) 

ID Cy^ O 



F2 = 9^ivv9^2vv-r^ FB{Mz,0,Mz,Mz,s,t) 
F^ = 9^,AV9^2AV ^ Z 4 FB{Mz,0,Mz,MA,s,t) 



16 c' 



w 



1 



+ 9^iH±W± 9^2H±W± -^Fb{Mw:0, Mw, Mh±, s,t) 
9^iAv 9^2AV 7^ Fb{Mz, 0, Mz, Ma, s, t) 



(A2) 
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where the couphngs g-i>vV} 9<!>az and gq>H±w are given in Table la. In terms of the Dijk 
four-point Passarino-Veltman functions [0], the form factors Fa^b = -^a.bI"^!, '"^2, f^s, ^i, 
s,t) [the analogous of Fi and F2 of eqs.(7) and (8)] are given by 



Fa = Dn + {u^t , M^, ^ M$J 



(A3) 



2(M4 -u)Do + 2(M|^ - u) Dn + 2(M4 -t)D,2 + 2t D13 (A4) 
1 D23 + 2sD2A + 2{t-Ml) D25 + \ ^ D2e + 8 D27 



t- M. 



D37 + 



u- M. 



2 2 

+ {u^t , ^ M$J 



with 



F>ijk = Ajfc("^l,"^2,^^3,^^4,0,0, M$^,M$2,s,t) 

For the class of diagrams consisting of sneutrino/chargino loops, the contributing 
amplitudes F^^ g are given by 



n 



^ Ml 

" 2^ '^~T- dex+y 9ex+o 9^iuu g^2vu FA{m;^, -m^+,mo, m^^, s, t) 
i=l ° 



■ 9ey+u 9e 



,^=^1 2CH."^^'^"^^^-^ 
J2 [d^^i^^ dl^xtx- ^(7"(-"^x+'"^^''^^'~"^X+'^'^'^*i'^*2' 



a=L,R 



and 



^ ^exfi' 9exti^ 



i,j,k=l 



a,b=L,R 



-^4,5.6 - 



(A6) 



Fa is given in eq.(A3), and the new form factors F^°' = F^"(mi, m2, ms, 7714, s, t, M$j, M^j) 
and = F^"*(mi, m2, ms, 7714, s, t, M$j, M^J, with a,b = L, R, are given by 



-R 



mi Dq + rriiD 



12 
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rp-RR 



■i^-RL 
D 



i^-LR 
D 



-LL 



TT+LR 



TT+RL 



p+RR 



2 "^4 Di2 

s (D12 + D13) + {s- MlJ D23 + s {D24 - D25 - D26) + 2 D27 
Ml^ L'as -{t- MlJ D37 -{u- Ml^) D39 - s ^310 - 6 L'313^ 



- mi m3 ( i:)o + -D13 ) 

1713 1714 Di3 



(A7) 



Finally, for the class of diagrams consisting of selectron/neutralino loops, the ampli- 
tudes -^7^8,9 given by 



Ml 



X/ 2 9ex'ieR 9ex'ieR 9'S>ieReR 9<S>2eReR FA{rneji,rn^o,meji.,meii, S,t) 
i=l ' ' 



and 



>^ Mz 

1^ J— 9ex°eR 9ex°~eR 



[9^ieReR 9l^x°X° ^c"("^x°' ^^R^ "^X°' ^' ^' ^$2) 

a=L,R 

+ fl'*2eHefl 5'$ix°x° Fc°'{m^o,mei„ mil,, "^x?> ^> -^*2> 

^ 1 

23 O ^ex»eH i/exOeH 



i,j,k=l 



X] [5'*ix«xO 5'4-2x;^x';! FD°'^{m^o,mej,, m^o,m^o, s, t, M$^, M^J 

a,b=L,il 

+ 9l2x°x°j 9iix°xl ^D°'''{iT^xp '^eR^m^o^rn^o.: u, M^^, M$j] (A8) 



-^7,8,9 - Ft^sA^R , F^,D F^j)) 



(A9) 



where the form factors Fa,b,c,d have been given previously. The normalized couplings of 
the neutral Higgs bosons to charginos, neutralinos and sleptons are displayed in Tab. lb 
and Ic. The only remaining couplings which have to be defined, are the electron- 
chargino-sneutrino and electron-neutralino-selectron couplings, which, normalized to 

r ^ -|l/2 

9w — v'^'Gf Mw, are given by 

9ex+i' = ' dex'^CR = "2 — iVjs , Qex^ei^ = -Ni2 - —N^ (AlO) 
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where the matrices V and N can be found in Ref.[^ 



For the pair production of the CP-odd Higgs boson in the MSSM, e^e — ^ AA, the 
differential cross section takes the form 



da 



dcos6' 



16 S TT^ 

(1 + A_)(1-A 



1 - 4^ {ut - M\) 



12 



fc=10 



+ (1-A_)(1 + A_ 



12 
k=10 



(All) 



where we have used the same notation as previously. The three contributing amplitudes 
-^10,11,12 are given by 
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glAyFB{Mz,^,Mz,Mh,s,t) 

1 



+ gHAvFB{Mz,Q,Mz,MH,s,t) + 



10 



gl^yFsiMz, 0, Mz, M;,, s, t) + gl^yPsiMz, 0, M^, M^, s, t) 



-^it — E 9 9ex°eR 9exleR 



11 



i,j,k=l 



E [^Ax«xO ^Ax«x^ F^''^{m^o, mg^, m^o, m^o, s, t, M^, M^) 



a,b=L,R 



+ ^Ax«xO ^Ax^x^ F^''\m^o,me^, m^o, m^o, s, m, Ma, Ma) 



^ 1 

E 2 9ex°eL 9exleL 
i,j,k=l 



a,b=L,R 



+ gAx1x° 3ax°xI Fd'^^^x^^'^^l^ '^X^ "^X°' 



-^12 



12 



i,j,k=l 



a,b=L,R 



(A12) 

with Fb and Fd as given previously; the couplings of the pseudoscalar Higgs bosons to 
the neutralino and chargino states are again given in Table Ic. 
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Table 1: 





g^vv 


g>i>AZ 




h 


sin(/5 — a) 


cos(/3 — a) 


=F cos(/3 — a) 


H 


cos(/3 — a) 


— sin(/3 — a) 


± sin(/3 — a) 


A 








1 



Tab. la: The Higgs-vector boson couplings g^vv [normalized the SM Higgs coupling Qh^vv = 
r 1 1/2 

2 \plGp -^y]) and the Higgs-Higgs-vector boson couplings [normalized to = (v^Gi;')^/^ 

Mvi/ and gz = {\/2GfY^'^ Mz for the charged/neutral weak couplings]; the latter come with the 
sum of the Higgs momenta entering and leaving the vertices. 



lilj 


ghijj 


gmjj 


gAijj 


eL^L 


(2s^ - 1) sin(/3 + a) 


-{2sl, - l)cos(/3 + a) 





eRCR 


sm{(3 + a) 


— 2s^ cos(/? + a) 





i>Li>L 


sm{P + a) 


— cos(/9 + a) 






Tab. lb: The couplings of the neutral Higgs bosons to left- and right-handed sleptons, nor- 

,1/2 



malized to g\y = V2Gf 



gL,R j ^ 


h 


H 


A 


„R 

^^xtx~ 


Q*i sin a — S*^ cos a 
Qij sin a — Sij cos a 


—Q*ji cos a — S*i sin a 
—Qij cos a — Sij sin a 


-Q*i sin f3 - S*i cos (3 
Qij sin P + Sij cos /5 


— 7^ — — 

3^x1x° 


Q'j* sin a + Sj* cos a 
Qij sin a + S'-j cos a 


—Q"ji cos a + Sj* sin a 
—Q"ij cos a + S'-j sin a 


-Q'j* sin (3 + Sj* cos /5 
sin /5 — S'-j cos /? 



Tab. Ic: The couplings of the neutral Higgs bosons to charginos and neutralinos, normalized 

r ^ -|l/2 „ „ ^ 

to gw = v2Gf the matrix elements Qij /Sij and Q^j/S^j can be found in Ref. [|21j] . 
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Figure Captions 



Fig. 1: Feynman diagrams contributing to the Higgs boson pair production process in e^e 
collisions: (a) SM Higgs production, CP-even (b) and CP-odd (c) Higgs boson 
production in the MSSM. 

Fig. 2: a) The cross sections for the Higgs pair production in the SM, e'^e~ — > H^H'^, 
as a function of the Higgs mass for two center-of-mass energies, ^/s = 500 GeV 
[dashed hues] and ^/s = 1.5 TeV [solid hues]. The lower curves correspond to 
the unpolarized cross sections, the middle curves to the cross sections where the 
electron beam is longitudinally polarized and the upper curves to the cross sections 
where both the electron and positron beams arc longitudinally polarized, (b) The 
angular distribution as a function of the scattering angle in the unpolarized case, at 
y/s = 500 GeV and for Mh = 150 GeV. 

Fig. 3: The cross sections for the pair production of the lightest Higgs boson in the MSSM, 
e''"e~ — > /i/i, as a function of for two center-of-mass energies, ^/s = 500 GeV and 
a/s — 1.5 TeV and for two values of tan (3 — 1.5 and 50. The solid curves correspond 
to the full cross sections, while the dashed curves correspond to the cross sections 
without the SUSY contributions. 

Fig. 4: The cross sections for the process e+e" HH, as a function of Mh for two center- 
of-mass energies, ^/s — 500 GeV and ^/s = 1.5 TeV and for two values of tan /3 = 1.5 
and 50. The solid curves correspond to the full cross sections, while the dashed 
curves correspond to the cross sections without the SUSY contributions. 

Fig. 5: The cross sections for the process e+e" AA, as a function of Ma for two center-of- 
mass energies, ^/s = 500 GeV and ^/s = 1.5 TeV and for two values of tan (3 = 1.5 
and 50. The solid curves correspond to the full cross sections, while the dashed 
curves correspond to the cross sections without the SUSY contributions. 

Fig. 6: The cross sections for the process e+e~ Hh, as a function of Mh for ^/s = 500 
GeV and for two values of tan/9 = 1.5 and 50. The solid curves correspond to the 
full cross sections, while the dashed curves correspond to the cross sections without 
the SUSY contributions. 
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